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Abstract. A general theory of thermal magnetic fluctuations near the 
surface of conductive and/or magnetically permeable slabs is developed; 
such fluctuations are the magnetic analog of Johnson voltage noise. 
Starting with the fluctuation-dissipation theorem and Maxwell's equa- 
tions, a closed-form expression for the magnetic noise spectral density is 
derived. Quantum decoherence, as induced by thermal magnetic noise, is 
analyzed via the independent oscillator heat bath model of Ford, Lewis, 
and O'Connell. The resulting quantum Langevin equations yield closed- 
form expressions for the spin relaxation times T\, T2, and Tip. For real- 
istic experiments in atomic physics, quantum computing, and magnetic 
resonance force microcopy (MRFM), the predicted relaxation rates are 
rapid enough that substantial experimental care must be taken to mini- 
mize them. At zero temperature, the quantum entanglement between a 
spin state and a thermal reservoir is computed. The same Hamiltonian 
matrix elements that govern fluctuation and dissipation are shown to 
also govern entanglement and renormalization, and a specific example 
of a fluctuation-dissipation-entanglement theorem is constructed. We 
postulate that this theorem is independent of the detailed structure of 
thermal reservoirs, and therefore expresses a general thermodynamic 
principle. 
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1. Introduction 

The engineering control of quantum decoherence is central to several emerg- 
ing scientific goals. One such goal is the direct observation of molecular 
structure — in three dimensions and with single-atom resolution — by mag- 



netic resonance force microscopy |25|. Another is the solution of otherwise 



intractable mathematical problems — like factoring large numbers — by quan- 
tum computing [p^] . 

At present, the main practical obstacle to achieving these goals is that 
the various mechanisms that cause quantum decoherence are not yet fully 
cataloged and understood. In consequence, quantum-coherent experiments 
often reveal unexpectedly fast decoherence rates, due to unanticipated — or 
even previously unknown — relaxation mechanisms. Similar struggles with 
decoherence have occupied previous scientific generations; Orbach for 
example reviews the forty-year struggle to achieve a reasonably comprehen- 
sive understanding of relaxation in electron spin resonance. 

This article is concerned with a decoherence mechanism which is ubiqui- 
tous in quantum-coherent technologies: thermal magnetic noise. Physically 
speaking, thermal magnetic noise is created by the same thermal current 
fluctuations that create Johnson noise. Thus, thermal magnetic noise is 
present in any device that contains electric conductors. 

We will be mainly concerned with practical engineering aspects of ther- 
mal magnetic noise, but we will also give an explicit example — the first in 
the literature to our knowledge — of a fluctuation-dissipation-entanglement 
theorem. Such theorems may be broadly defined as invertible functional re- 
lations between the dissipative kernel of a system and the system's quantum 
entanglement with a thermal reservoir, such that the entanglement deter- 
mines the dissipation and vice versa. 

This article is organized as follows. Prior work relating to thermal mag- 



netic noise is reviewed in Section 1.1. The main new results are summa- 



rized in Section |2[ these results consist of closed-form expressions for the 
magnetic noise spectral density (Section p.l[ ), quantum decoherence (Sec- 
tion |2.2[ ), and a fluctuation-dissipation-entanglement theorem (Section |2.3D . 
The practical implications of these results for quantum-coherent engineering 
are discussed in Section ^ via worked examples for trapped atom experiments 
(Section |3.1D , magnetic resonance force microscopy (Section 3^), and quan- 



tum computing (Section |3.3D . Some topics for further research are suggested 
in Section ^. 

Algebraically tedious derivations are relegated to three appendices. Ap- 
pendix ^ derives a closed-form expression for the magnetic noise spectral 
density. Appendix |B| derives quantum Langevin and Bloch equations that 
describe decoherence in the presence of thermal magnetic noise, and Appen- 
dix ^derives fluctuation-dissipation-entanglement theorems for both spin-^ 
particles and harmonic oscillators, as well as relations between dissipation 
and renormalization. 
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1.1. Prior work relating to thermal magnetic noise. The quantum 
field theory of electromagnetic fluctuations in lossy linear media is presented 
in textbooks by Landau, Lifshitz, and Piaevskii |l^ and by Rytov et al. p^ . 
Far from a warm body, this theory describes fluctuations which are simply 
the familiar phenomenon of black-body radiation. Near to a warm body, 
the situation is considerably more complex. Non-radiating terms in the field 
theory become important, and create phenomena such as the Casimir force 
(identical to the Van der Waals force), which is familiar from chemistry and 
force microscopy. Recently, Dorofeyev et al. Q have observed the attractive, 
fluctuating and dissipative components of the Casimir force, which were in 
good agreement with field-theoretic predictions. 

The Casimir force arises mainly from fluctuating electric fields, which do 
not directly couple to the spin magnetic moments. Our investigation instead 
focuses on thermal fluctuations in magnetic fields. 

To the best of our knowledge, the theory and experiments of Varpula 
and Poutanen as reviewed and extended by Nenonen, Montonen, and 



Katila [19|, were the first experimental observation and theoretical analy- 
sis of thermal magnetic fluctuations. Their work focussed on biomagnetism 
experiments, in which magnetic fluctuations originate in the copper walls 
of shielded rooms; such noise can be easily large enough to obscure bio- 
magnetic signals. They presented a phenomenological model of thermal 
magnetic noise in which metallic conductors were modeled as collections of 
independent thermally excited resistive elements. The resulting predictions 
were in excellent accord with experiment. 

Quantum-coherent technologies operate in a vastly different regime from 
biomagnetism experiments. Ambient temperatures T are as low as is ex- 
perimentally feasible — a few Kelvins down to millikelvins — and observation 
frequencies oj are megahertz to gigahertz, such that ^ k^T. 

Quantum-coherent phenomena dominate this physical regime. The ap- 
plicability of the Varpula-Poutanen noise model then becomes uncertain, 
because the model is derived from a purely phenomenological description of 
thermal noise in room-temperature conductors. We have therefore sought 
to improve and extend the Varpula-Poutanen model in five respects: 



1. Rigorous quantum mechanical and thermodynamic foundations have 
been provided via the fluctuation-dissipation theorem. 

2. The results now encompass superconducting materials and/or materi- 
als with non-vanishing magnetic susceptibility. 

3. A closed- form expression for the thermal magnetic noise spectral den- 
sity has been obtained. 

4. Bloch equations describing quantum decoherence have been derived 
via a quantum Langevin formalism and the independent oscillator heat 
bath model of Ford, Lewis, and O'Connell Q. 

5. Fluctuation-dissipation-entanglement theorems have been constructed, 
and the link between dissipation and renormalization has been clarified. 
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In carrying through this program, we have been able to prove that the 
Varpula-Poutanen model is exact in the physical regime relevant to biomag- 
netism: real conductivity, macroscopic length scales, and large temperature. 
This explains why the Varpula-Poutanen model yields predictions which are 



in excellent accord with experimental measurements of biomagnetism |26], 
and is a tribute to their pioneering physical insight. 

2. Physical motivation and main results 

The necessary existence of thermal magnetic noise can be deduced by 
considering a loop antenna held near a conducting slab, as shown in Fig. ^. 
If a voltmeter is placed across the terminals of such an antenna, a fluctuating 
voltage V{t) will be observed. The spectral density^ ^^/(a;) of this Johnson 
noise is related to the complex antenna impedance Z{lij) by the well-known 



equation |12 



(2.1) Svi^^) = Re[Z{uj)) fkj coth (^ 2^t } ' general; 

~ Re{Z{uj)) 2kBT, for kj^T » hw. 

Here Re(Z'(a;)) is the resistive impedance. This result is just the fluctuation- 
dissipation theorem as it applies to fluctuating voltages in dissipative imped- 
ances. 

Now we ask, what is the physical origin of these voltage fluctuations? 
We suppose that the coil has negligible intrinsic resistance, in which case 
the voltage can only have been induced by a magnetic flux threading the 
antenna loop. In turn, this magnetic flux can only have been generated by 
thermally excited currents in the nearby conducting slab. 

So wherever Johnson noise is present, thermal magnetic noise will be 
present also, because the same underlying current fluctuations generate both 
phenomena. 

2.1. Thermal magnetic spectral densities. We will begin by presenting 
closed-form expressions for the spectral density of thermal magnetic noise, 
and discussing their physical significance. Deriving these expressions is te- 
dious, and is deferred until Appendix |A[ 

Since the thermal magnetic field B{t) is a vector, its autocorrelation 
Cb{t) and spectral density Ssi'^) are matrix-valued functions defined by 

(2.2) CB(r) = (S(0)®B(r)), 

/•oo 

(2.3) 5b(c^)= / e~'^^CB{T)dT 



^Our normalization for spectral densities is Sv{lo) ~ f^^dre ^'^^ {V{0)V{t)), where 
{V{0)V{t)) is the voltage autocorrelation. Thus {V'^) — l/(27r) J^^duj Sv{i^), so that 
spectral densities are two-sided (positive and negative frequencies) , with bandwidths given 
in Hertz. 
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where O denotes the outer product of two vectors: (a b)ij = aibj. As 



shown in Section A^, S'b(<-^) is of the general form: 



(2.4) 5BH=rH?^coth(^^) 

= ["/ + n (g) nl r(u;);iu;coth ( f ^ ] 

\2kBT J 

where h is the unit vector normal to the slab surface and / is the 3x3 
identity matrix. 

The magnetic dissipation tensor T{uj) in (|2.4| ) plays the same fundamen- 
tal role that resistance plays in Johnson noise, in the sense that — as we 
will see — knowledge of r(a;) quantitatively determines a broad spectrum of 
fluctuation, dissipation, and quantum entanglement phenomena. 



2.1.1. Device parameters. In general, we will concerned with thermal mag- 
netic noise induced at a distance d from a conductive slab of thickness t, 
as illustrated in Fig. |^. The slab is assumed to have a frequency-dependent 
complex conductivity cr{u;) and permeability fJ-{uj); for compactness we sup- 
press the frequency-dependence of a and fi. All results are presented in S.I. 
units, and all time-dependence is assumed to be e*"^*. See Section [A.3| for a 
statement of Maxwell's equations in this convention. 

The four device parameters {d,t,a,fi} may be regarded as fundamental; 
all other parameters are derived from them. The single most important 
derived parameter is the skin depth A = Icvfial"^^"^ . Another important 
derived parameter is the phase (j) of the complex conductivity, as defined by 
j = aE = \a\e~'^'f', where E is the E- field and j is the current. 

2.1.2. Superconductors. For the special case of an ideal superconductor, the 
London equation can be stated in the form i^QUjj = A^-E, where Al is the 
London penetration depth. Comparing this with the definition of the con- 
ductivity j = aE, we see that for an ideal superconductor a = —iXl/{fioUj), 
A = Al, and (j) = tt/2. Thus ideal superconductors do not have infinite 
conductivity, but rather have a wholly imaginary conductivity. 

Real superconductors approximate ideal superconductors only at zero 
temperature. At finite temperatures a mixture of superconducting and nor- 
mal phase charge carriers is present, such that the conductivity has an inter- 
mediate phase (p G (0)7r/2). See Q for a review of the literature, including 
measurements of a at finite temperature. 

2.1.3. Single paramagnetic or diamagnetic conducting slabs. Most nonfer- 
rous metals and insulators have relative magnetic permeability K = /////q — 
1, where // is the permeability of the material (the slab in our case), and //q 
is the permeability of the vacuum. In Appendix ^ we derive the following 
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asymptotic dependance of the magnetic dissipation for if = 1 

t 



,2.5) V(.) = >i^.l^ 



d{t + d) ' 
3A3 



d4 cos(7r/4 - 0/2) 



A » min{d, v dt}\ 
A <C min{d, t}; 



2A^(dt-8A2sin</.) , 
^5^2 ^- 



Collectively, these three limits span the entire {d,t^a,uj} design space. The 
following expression smoothly interpolates all three limits: 

(2 61 Tiuj) ~ 3^gRe(o-)A^tcos0 

^ 647r(3A3(i((i + t) + (l-e-"0td2(^ + 2A)2cos(7r/4-(/./2))' 



It is constructed by adding the first and second asymptotic limits of ( |2.5| ) 
inversely, then introducing a transition parameter 

(2.7) Oc = ((it + 8A2sin0)/(2dAcos(7r/4-(/)/2)) 

such that the third asymptotic limit also is smoothly interpolated. 

A numerical survey shows that the resulting closed-form expression ( |2.6| ) 
predicts noise levels that are no more than 1.75 dB too large and 0.5 dB too 
small everywhere in {d, t, <t, a;} design space — an accuracy which suffices for 
almost all practical design work. 

2.1.4. The case of two conducting slabs. The case of magnetic fluctuations 
measured at a point midway between two identical conducting slabs is also 
solved in Appendix|A|, and the resulting asymptotic expressions are precisely 
double those given in ( p. 5]) . A numerical survey shows that doubling the one- 
slab noise yields a slight overestimate of the two-slab midpoint noise, but 
the error is no more than 0.6 dB for all {d,t,a,uj}. Thus, to a very good 
approximation two adjacent slabs each generate independent noise at the 
midpoint. 

2.1.5. The general case. In the general case in which both lossy magnetic 
permeability and electric conductivity are combined with a skin depth A 
comparable to d, no analytic results are presently available. Appendix ^ 
presents a one-dimensional integral (A. 18) which must be numerically eval- 



uated to determine the magnetic spectral densities. 

2.1.6. Comparison with experiment. Although our general integral expres- 
sion ( |A.18| ) is superficially of a different form than the Varpula-Poutanen 
model [^6[, and is derived from a very different underlying model, a numer- 
ical comparison shows that these two models yield identical predictions in 
the physical regime relevant to biomagnetic experiments: real conductivity, 
audio frequencies, vanishing magnetic susceptibility, and large temperature. 
Since in this regime the Varpula-Poutanen model has been shown to agree 
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very well with experiment, we may regard the model presented in this article 
as having been at least partially validated. 

More extensive comparisons of our model with the Varpula-Poutanen 
model are not feasible, because the V-P model does not contain parame- 
ters corresponding to the magnetic susceptibility fi, conductivity phase (p, 
or the quantum of action h. 

We note also that no model of thermal magnetic noise has yet been ex- 
perimentally tested at microscopic scales, low temperatures, and high fre- 
quencies, in materials that are wholly or partially superconducting, or are 
magnetic — which is precisely the regime of greatest importance to quantum- 
coherent engineering. 

2.2. Quantum decoherence. Now we turn our attention to the quantum 
decoherence induced by thermal magnetic fields. We assume a two-state 
system coupled to the magnetic thermal field B{t) via the Hamiltonian 

H = --f{Bo + B{t))-s 

where 7 is the gyromagnetic ratio of the spin, Bq is a constant polarizing 
field, and the spin angular momentum s satisfies the usual commutation 
relations [si, Sj] = iheijk Sfc. 

Our results will be expressed in terms of the rotating-frame magnetic 
spectral density S^Brot('^)> which is given in terms of the laboratory-frame 
density Ssi^) ( |2.4| ) via 

(2.8) Sb^M = {b®b) tT[b^b- Sb{^)] + 

(J- 60 6) ^tr[(/-6®6).SB(u;o)], 

where the precession frequency ujq about the unit axis of spin precession 6 
is given by 7.60 = ^^ob, and tr[. . . ] is the matrix trace (the sum of diagonal 
elements) . 

In Appendix ^ we derive the quantum Langevin equations for a spin-^ 
particle coupled to a thermal magnetic field, and from them we show that 
the time evolution of the spin density matrix is described by Bloch-type 
equations with relaxation times Ti, T2, and Tip given by 

(2.9a) 1 = i 7' tr [(/ - 6 ® 6) . (0)] 

(2.9b) 1 = _L + [6 ® 6 . 5b.., (0)] 

(2.9c) _L = i^2t^[(^ _ ^ . Sb^M] 

Here Tip is the relaxation time in the presence of a radio-frequency (RF) ap- 
plied field — commonly called a "spin-locking" field. In the frame co-rotating 
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with the RF field, it appears as a constant B-field Bi = Bibi whose char- 
acteristic precession frequency is uji = 7^1. Substituting dU-U) in ([2!9al - 



2.9b ), we obtain {Ti,T2,Tip} exphcitly in terms of the magnetic damping 



coefficient r(a;) given in (|2.6|) 



(2. 10a) 



— - 7^ (3 — cos^ 9) T(ujr\) hujQ coth ( — — ^ | 



(2.10b) l- = ±. + \^\i + cos'9)T{Q)2kBT 



1 _ 1 + cos2 (3 

sin^P (1 + cos2 9) T{u,) fu^i coth 

Here cos 9 = b • h and cos /? = b • bi, where bi is the RF field axis in the 
rotating frame, b is the polarization axis, and h is the unit vector normal 
to the slab surface. 

Note that r(a;) appears with varying frequency arguments lo £ {0, wi, wq}; 
this has important design consequences because in most cases T{uj) exhibits 
a strong frequency dependence. 

As far as practical quantum engineering is concerned, expressions (|2.10a| - 



2.10c ) are the main "deliverable" of this article. 



2.3. A fiuctuation-dissipation-entanglement theorem. From a purely 
formal point of view, fluctuation-dissipation-entanglement theorems exist for 
a simple reason: the same Hamiltonian matrix elements that control fluctua- 
tion and dissipation also control certain measures of quantum entanglement; 
Appendix ^ discusses this point of view. 

Physically speaking, we reason as follows. We consider a two-state quan- 
tum system, which — as in the real world — interacts weakly with a thermal 
reservoir. We adjust the temperature of the reservoir to zero, and ask: what 
is the probability that the two-state system is not in its ground state? To the 
extent that this probability is non-zero, it describes an irreducible quantum 
entanglement with the thermal reservoir. 

Now, at zero temperature the classical entanglement probability is zero. 



and even the Bloch equations ( B.19 ), which have a firmer quantum justi- 



fication, relax spin systems to their ground state at zero temperature, and 
thus predict zero entanglement. 

However, a higher-order calculation reveals that the entanglement prob- 
ability is finite even at zero temperature. We calculate this probability as 
follows. As before, we consider a spin-^ particle magnetically coupled to 
a thermal reservoir. We polarize the spin via an external field Bq, which 
lifts the degeneracy of the spin system along a polarization axis ujqP = 'jBq, 
with 7 the gyromagnetic ratio of the spin and (Jq the (by convention posi- 
tive) precession frequency. We let \ipo) be the ground state of the isolated 
spin, and we define a projection operator Pg = I — jV'o) ® {tpo\, with I the 
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identity operator. By construction, the expectation value of Pg vanishes for 
an isolated ground state: {iPo\Pe\iPo) = 0. 

Of course, no real-world spin system is perfectly isolated.^ We therefore 
enlarge our Hilbert space to encompass a thermal reservoir with a potentially 
infinite set of basis states i G 0, 1, . . . , oo}. We generalize to include 
the thermal reservoir by defining a new projection operator £ 

(2.11) s. = ^\<l)i)<S)PE0{<pi\ 

i 

which projects onto the subspace in which the spin is not in its ground state. 

It is a well-defined problem to calculate, perturbatively, the expectation 
(^'o|£|^o)) where |^'o) the ground state of the combined spin-plus-thermal- 
reservoir system. Keeping in mind that £ and |^'o) depend on p and ojq, 
we can define a scalar entanglement function E{p,ujq) = (^'o|£|^'o)- In 
Appendix we show that E{p,ujq) is related to the magnetic dissipation 
tensor T{uj) by the following fluctuation- dissipation- entanglement theorem: 

(2.12) E{p,uJo) = ^ 02{c^tr[(/-p(g)p) .r(w)];a;o} 

Here ©2 is a Stieltjies transform, in Bateman's notation ^j, as defined by 

/•oo 

(2.13) 6p{/(x);y}= / dx f{x){x + y)~P 

Jo 

We note that Stieltjies transforms are invertible, such that a measurement 
of the entanglement function E{p,ujq) suffices in principle to determine the 
magnetic dissipation tensor r(a;) and vice versa. 

For practical quantum engineering purposes, the main utility of this the- 
orem is that it allows us to determine the hard-to-measure entanglement 
from the easier-to-measure (or predict) dissipation. For the particular case 
of thermal magnetic noise, the frequency dependence ( ^.6| ) of r(a;) is such 
that the Stieltjies integral is absolutely convergent, so that the entanglement 
of real-world devices can be readily be calculated. 

2.3.1. Approximate expression for the entanglement. In general, the Stieltj ies 
integral ( |2.12| ) cannot be evaluated in closed form. However, a simple and 
physically meaningful approximate evaluation is possible. We begin by re- 
marking that dissipative kernels always have a cut-off frequency Wc, because 
otherwise their associated noise spectral density would carry infinite power. 
For the particular case of thermal magnetic noise originating in a thick plate, 
the cut-off frequency is such that A(a;c) ~ d. We assume the frequency of 
interest loq is small compared to the cut-off frequency, so that ^ wq, and 
we brutally approximate r(a;) as constant over the range uj G (wc^c); and 

^ If nothing else, the spin can exchange quanta with the vacuum, which may be re- 
garded as a zero-temperature thermal reservoir. It follows that fluctuation-dissipation- 
entanglement theorems can be constructed in quantum field theories, where their signifi- 
cance remains to be elucidated. 
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zero elsewhere. Then the Stieltjies integral ( |2.12| ) can be evaluated in closed 

form as 

(2.14) 

E{p, ujo) ~ i— tr [(/ - p ® p) • r(wo)] InKM) from 



4tt 

f^tT[{I-p0P). Hlu./uo) by (U for T ^ 



^""^^^/^^^ by (U) and (S) 



27ru;o^i 

where Ti is evaluated at zero temperature. 



3. Quantum engineering implications 



We will now apply (^.4|- 2.10| ) in the design analysis of representative quan- 



tum technologies. Our goal is partly to illustrate practical applications of 
our results, and partly to identify areas where further research is needed. 

A small but important point: when quoting numerical values oc 1 /\/Hz we 
shall embrace the usual engineering convention that bandwidths encompass 
only positive frequencies; this requires the insertion of an additional factor 
of two when evaluating two-sided spectral densities, e.g. (|2.4|). 



3.1. Atom and ion traps. We begin by considering thermal magnetic 
noise at centimeter length scales and audio frequencies, at room tempera- 



ture. This is the same regime considered by Varpula and Poutanen |19,26|, 
and it is also a regime typical of at least some atomic physics experiments [|]. 

Specifically, we will calculate the spectral density of the thermal magnetic 
fields between room-temperature two copper slabs, each t = 1 cm thick, 
spaced 2d = 2 cm apart. We are particularly interested in the zero- frequency 
spectral density in the n direction normal to the slab surface, because this 
describes thermally induced fluctuations in the background polarizing field. 
At room temperature T = 300 K the conductivity of high-purity copper 
is of order 5.9 x lO''' m)~^ From ( |2^ ) and (|2^), with d = t = 1 cm 
and (j) = 0, and taking into account that we have two plates whose noise 
is additive, we find at zero frequency (n • 5^(0 Hz) • n)^/^ ~ 1.2 pT/a/Hz, 
rolling off to (n • Sb(IOOHz) • n)^/^ ^ q g pX/^H^ at 100 Hertz. 

Until recently, such picoTesla fluctuations would have been viewed as be- 
ing of no practical consequence. However, magnetic fields change sign under 
time reversal, and hence magnetic fluctuations locally violate time-reversal 
invariance, and so these fluctuations must be understood and controlled in 
high-precision tests of fundamental physics. 

For example, the most stringent experimental limit on the electric dipole 
moment d^ of a fundamental particle is d^ < 8.7x 10"^*^ e cm for the mercury 
isotope ^^^Hg Q (here e is the electron charge). Typically, dipole moments 
are measured by observing the change in precession frequency induced by an 
applied electric field of magnitude E. It follows that S b can be expressed 
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as an equivalent dipole noise spectral density Sd = {hj/ h • Sb{0) • n. 
For ^^^Hg, the gyromagnetic ratio 7/(27r) = 7.59 MHz/Tesla. This yields a 
zero-frequency equivalent dipole noise for our example of room-temperature 
copper plates of 

ES]''^ ~ 1.94 X 10^20 eV/^. 

For a typical electric field E ~ 10'^ V/cm, the equivalent dipole noise would 
therefore be 5"^^ ~ 1.94 x 10^24 g cm/^/lfe. 

This is a substantial noise level: it would naively require an averaging time 
on the order of 10^ seconds to better the published electric dipole moment 
limit . And it would not be entirely straightforward to reduce the thermal 
magnetic noise by making the copper shields thinner or less conductive, 
because this would defeat their purpose of shielding the experiment from 
external fields. 

Fortunately, another effect provides mitigation: trapped atom experi- 
ments typically measure the net signal from atoms in many different regions 
of a cell. To the extent that n independent regions are averaged, the equiv- 
alent dipole noise power will be reduced by a factor of 1/n. Some of the 



formalism necessary for calculating n is set forth in Appendix A.l. Detailed 



calculations would be strongly dependent on the particular design chosen. 



3.2. Magnetic resonance force microscopy. Magnetic resonance force 
microscopy (MRFM) is a quantum-coherent technology whose ultimate ob- 
jective is to produce magnetic resonance images of individual molecules in 
situ, nondestructively, in three dimensions, with Angstrom resolution. Such 
a technology would allow much of molecular biology to be conducted as an 
observational science — along the lines of, e.g., astronomy — rather than an 
experimental science. 

We remark that a single-spin MRFM imaging device can be alternatively 
regarded as a first-generation solid-state quantum computer, in which the 
individual qubits carry binary information about the presence or absence of 
a spin spin at a specified atomic coordinate. 

There are two main experimental challenges in MRFM. The first is that 
single-spin signal forces are exceedingly small, of order 10~^^ N for electron 
magnetic moments. The MRFM community's design options for achieving 
the required sensitivity are reasonably well understood (but challenging to 
implement in practice): reduce the mass of the cantilever, increase its damp- 
ing time, reduce the ambient temperature, and employ a sharper magnetic 
tip with a higher field gradient. 

The second — emerging — challenge in MRFM is that the spin state must 
maintain its quantum coherence for a time long enough to detect the signal 
force. Here the design issues are not yet well understood. We will illustrate 
some of these issues by calculating the effects of thermal magnetic noise on 
spin relaxation in a typical MRFM environment. 
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We model the magnetic tip as a sphere with radius r = 1 and uniform 
magnetization ^qM = 1 T. The electron magnetic moment is located at a 
separation d = 50 nm from the tip surface. At this distance the polarizing 
field is Bq = 2fiQMr^/ (3(r+(i)^) = 0.58 T, and the spin precession frequency 
is jBo/{2tt) = 16.1 GHz, where j/{2tt) = 2.8 MHz/T is the gyromagnetic 
ration of the electron. We take the ambient temperature to be T = 4 K. 
We note that d <^r and we therefore — roughly — model the tip as a slab of 
thickness t = d. Finally, the direction of the polarizing field is taken to be 
parallel to the vector normal to the tip suface {9 = 0). 

With the parameters {d,t,T,j,uJo,0} now specified (see Section 2.1.1 for 
further discussion), the electron spin relaxation rates 1/Ti{a) and 1/T2{a) 
are predicted as functions of the tip conductivity a by expressions (|2.10a| - 
2.10b| ), with results as shown in Fig. ^. For completeness, relaxation times 
for proton magnetic moments at the same distance from the tip are also 
shown. 



3.2.1. Design lessons for MRFM. It is clear that tip conductivity is a key 
engineering parameter in MRFM. For example, if the tip conductivity is 
as great as that of pure elemental iron (~ 4 x 10^ ri-m at 4 K) or nickel 
(~ 1.7 X 10^° f^-m at 4 K), the predicted relaxation rate is 1/Ti > 10^ s^^. 
Such rates are much too rapid to permit coherent spin imaging. 

Real- world MRFM tips typically are composed not of pure elements, but 
of magnetically "hard" alloys such as samarium-cobalt or neodymium-iron- 
boron. These alloys are commerically prepared as powders, without much 
regard for minimizing lattice defects: the resulting conductivity is very likely 
to be substantially reduced relative to the values quoted above for pure 
elements. There is not much data in the literature relating to the electrical 
conductivity of ferromagnetic particles at cryogenic temperature. Misra et 
al. [17] have measured the conductivity of sputtered Cu/Cr multilayered 



films to be ~ 2.5 x 10~^ VL-m in the zero-temperature limit. If we take this 
value to be representative of micron-scale ferromagnetic tips at cryogenic 
temperatures — in the absence of better data — the predicted relaxation rate 
is 1/Ti ~ 10 s~^. This is still an uncomfortably rapid relaxation rate. 

If these predictions are correct, the MRFM community has at least five 
design options: (1) reduce the temperature, (2) reduce the tip conductivity 
{e.g., by employing an insulating magnetic garnet tip), (3) employ super- 
conducting tips {e.g.., generating gradients via trapped vortices), (4) detect 
electron moments incoherently, or (5) focus on coherent proton detection, 
for which the predicted relaxation times are much longer. 

We see how vital it is to achieve a thorough understanding of all the 
mechanisms that influence spin relaxation. It is sobering, therefore, to re- 
flect upon some of the thermal reservoirs that we have not mentioned, which 
surely are present in ferromagnetic tips: spin waves, thermally excited do- 
main wall motions, and unpaired electron spins in passivating oxide layers, 
to name three. As discussed in Section ^, much work remains to be done. 
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3.3. Quantum computing. Kane has described a design for a solid- 
state quantum computer in which the qubits are the electron and nuclear 
spin quantum numbers of ^^P donor sites in a silicon lattice. Now we will 
consider the effects of thermal magnetic noise within Kane's device. 

In Kane's design, the ^^P donor sites carry both nuclear and electron 
quantum numbers, which are coupled by a Hamiltonian Hi^t of the form 

(3.1) -f^int = 7e-Bo • Se + JnBo ■ S,, + AH'^A ■ 

Here s^. and Sn are the electron and ^^P spin operators, 7e and 7n are the 
gyromagnetic ratios (in a convention where both ratios are positive), and 
the hyperfine coupling A has the value hA/{2TT) = 29 MHz. 

The computational action occurs between the two lowest energy levels 
of this system, and the presence of the hyperfine coupling modifies these 
lowest-energy levels in two respects. First, the energy separation of the two 
levels depends on A via an equation given by Kane as 

(3.2) hwo = -fnBo + 2A + 0{A^) 

Second, the two lowest energy levels of the system couple to external time- 
dependent B-fields Bextit) via an effective Hamiltonian Hes which a straight- 
forward perturbative calculation shows to be 

(3.3) = ^ B,xt{t) ■ K • a 

where cr = {ax,cry, az} is a vector of Pauli matrices acting on the two lowest- 
energy states, and K is a dimensionless matrix given by 

(3.4) K = 6®b+ + +0(^2) 

where & is a unit vector in the direction of the applied B-field. We ob- 
serve that the hyperfine coupling A acts to increase the coupling of the two 
lowest-energy states to thermal magnetic noise. This increased coupling is 
most conveniently taken into account by the replacement in ( |2.8D of the 
laboratory-frame spectral density Sb by an effective spectral density 

(3.5) Sb^S'§ = K-Sb-K 

Subsequent calculations of relaxation rates are unaltered. In Kane's design, 
the renormalized coupling K increases the predicted relaxation rates by a 
factor of order (1 + {2A) / {-fnBo)f ~ 7.2 . 

Per Kane's design [P, we specify a polarizing field Bq = 2T and a temper- 
ature T = 100 mK. We model Kane's A-gates and J-gates as metallic pads 
of thickness t = 5 nm at a distance d = W nm from the ^^P centers. The 
direction of the polarizing field is taken to be parallel to the vector normal 
to the pad surface {8 = 0). The electron-spin and nuclear-spin relaxation 
rates 1/Ti and I/T2 are then predicted as functions of the pad conductivity 
a by expressions (|2.10a - 2".10b| ), with results as shown in Fig. ^. 
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3.3.1. Design lessons for solid-state quantum computing. Kane's article op- 
timistically quotes the magnetic resonance literature to the effect 

at temperatures T 1.5 K the electron spin relaxation time is 
thousands of seconds, and [ . . . ] at millikelvin temperatures the 
phonon limited '^^P relaxation time is of the order of 10^^ seconds, 
making this system ideal for quantum computation. 

We see that thermal magnetic noise originating in the A- and J-gates is 
likely to substantially increase these relaxation rates. Fortunately, to the 
extent that calculation occurs only between the two lowest-energy states, 
the predicted rapid relaxation of the higher-energy electron states may not 
pose a problem. And it is gratifying that the predicted ^^P relaxation rates 
remain reasonably small even when the amplifying effects of the hyperfine 
coupling are taken into account. 

These predictions should be regarded very cautiously, however. As we 
discuss at greater length in the following section, even the notion of conduc- 
tivity becomes suspect at these length scales and temperatures. In our view, 
the main design lesson is that all potential noise sources must be treated 
seriously. 



4. Conclusions 

We have developed a general theory of thermal magnetic fluctuations 
near the surface of conductive and/or magnetically permeable slabs. We 
have applied the resulting closed-form expression for the magnetic spectral 
density (Section |2.1D to practical problems in atomic physics experiments, 
magnetic resonance force microscopy, and solid-state quantum computing 
(Sections |3.1[-^^). The theory predicts magnetic noise levels that are large 
enough to require substantial care to minimize their effects. 

More generally, our formalism indicates that if the dissipative kernel of 
a linearly damped system is known, then the consequences for fluctuation, 
dissipation, entanglement, and renormalization are fully determined, via 
fluctuation-dissipation theorems (eqs. ( 2!^ ) and (|C.S|)), fluctuation-dissipa- 



tion-entanglement theorems (eqs. (|2.12 ) and ( p.lCp ), and renormalization 



dissipation relations (eqs. ( |C.14D and ( C.15D ). 



We postulate that these fluctuation-dissipation-entanglement-renormaliz- 
ation relations are universal, in the sense that they apply to all thermal 
reservoirs, without regard for the internal structure of the reservoir, with 
the sole restriction that the coupling to the reservoir must be weak enough 
to be linearly dissipative. 

A major limitation of our article is that we derived our results in the 
context of a specific model of thermal reservoirs — the independent oscillator 
(10) model — and hence we have not proved universality. Ford, Lewis, and 
O'Connell's pioneering article on independent oscillator models proves 
universality in the restricted sense that these models are shown to describe 
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all forms of fluctuation and dissipation that are causal and linearly dis- 
sipative. An important milestone for further research would be to prove 
universality in an expanded context which included entanglement and renor- 
malization in addition to fluctuation and dissipation. 

There are also substantial reasons to expect that measurement processes 
might also be subsumed under a suitably expanded notion of universality. 
We reason as follows. Quantum-coherent technologies necessarily include 
some means for reading-out binary data; in quantum computing the data 
are qubits, in MRFM experiments the data are the presence or absense 
of a spin at a given atomic coordinate. In most cases measurement is a 
continuous interferometric process, achieved e.g. via weak interactions of 
a cantilever with laser-supplied photons From an engineering point of 
view, the measurement photons comprise a "good" reservoir — effectively 
at zero temperature — whose finely tuned interactions compete with many 
different "bad" thermal reservoirs that contribute only noise. The quantum 
system seeks to achieve equilibrium with the "good" and "bad" reservoirs 
simultaneously; the resulting quantum dynamics are traditionally described 
by master equations of the Fokker-Planck type. A unified formalism of this 
sort, treating measurement and noise processes on a common footing, would 
offer many practical advantages in quantum-coherent engineering. 

In carrying through these calculations, we have come to appreciate an 
emerging common ground — the technical challenge of preserving and ma- 
nipulating quantum coherence — that is shared by magnetic resonance force 
microscopy and quantum computation. Single-spin MRFM imaging devices 
(if they can be demonstrated) would constitute first-generation solid-state 
quantum computers, in which the individual qubits carried binary infor- 
mation about the presence or absence of a spin spin at a specified atomic 
coordinate. 

Arguably the most daunting challenge in developing quantum-coherent 
technologies is achieving simultaneous engineering control of a panoply of 
thermal reservoirs: 

• phonons and surface vibrational modes 

• ferro- and ferrimagnetic spin waves and domain walls 

• conduction band excitations in metals and semiconductors 

• paramagnetic and nuclear spins 

• mobile molecules and charges on surfaces 

• lattice defects and impurities 

• super-conducting condensates admixed with normal-phase conductors 

This list includes most of the excitations that are commonly studied by 
condensed matter physicists. We remark it is not enough to achieve good 
engineering control of some or even most of these noise sources; all must 



■^We remark that the detailed quantum dynamics of the interferometer-cantilever-spin 
system remains to be worked out. 
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be controlled simultaneously, as even one poorly-controlled noise source is 
enough to destroy quantum coherence. 

The dissipative kernels of these excitations are poorly understood at the 
ultra-low temperatures and mesoscopic length scales that are characteristic 
of quantum-coherent technologies. For example, even Maxwell's equations 
in conducting media — an exceedingly well-studied area of physics — will re- 
quire re-interpretation due to our still-limited understanding of mesoscopic 
conductivity |[7|,p^. 

In view of these difficulties, Landauer has modestly proposed |^4| that 
the following disclaimer be appended to all quantum computing proposals: 

This proposal, like all proposals for quantum computation, relies 
upon speculative technology, does not in its current form take 
into account all possible sources of noise, unreliability and man- 
ufacturing error, and probably will not work. 

Our article may be read as a response to Landauer's concerns: we have 
analyzed thermal magnetic noise as a "possible source of noise and unrelia- 
bility" in MRFM imaging and in Kane's proposal for a quantum computer. 
More ambitiously, we have done so in the context of a general formalism 
which might be applied across a broad range of thermal reservoirs. In do- 
ing so, we hope to have contributed to the emerging task of providing solid 
and well-organized foundations for the development of quantum-coherent 
technologies. 
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Appendix A. Fluctuation and dissipation 

Many readers will be acquainted with the fluctuation-dissipation theorem 
as it applies to voltage noise — known as Johnson noise — in resistive circuit 



elements |12]. The theorem can be briefly stated as follows: if a net charge 
q{t) = qe'-^^ is passed through a frequency-dependent impedance Z{uj), such 
that the current i{t) = dq(t)/dt = iujqe^'^^, and the power P{oj) dissipated 
in the impedance is P{oj) = Re(Z(a;)) g^w^, then the spectrum of thermal 
voltage fluctuations is given in terms of the dissipative impedance Re(Z) by 
(P). 

The general fluctuation-dissipation theorem allows us to deduce, by an 
exact analogy, that if a time-dependent magnetic moment m{t) = me*'^* 
is adjacent to a conductive slab, such that the power dissipated in the slab 
is P = m • T{uj) • ma;^, then the spectrum of the magnetic fluctuations is 
given in terms of the magnetic dissipation matrix V{u}) via (|2.4D . 

Discussions which justify the above line of reasoning are presented in 
many graduate- level textbooks |j^,|ll[ , but these discussions tend to be rather 
abstract and lengthy. Some readers will prefer the following physical deriva- 
tion, which is completely general and rigorous. 

We will deduce the spectrum of magnetic fluctuations from the well-known 
spectrum of thermal voltage fluctuations (Johnson noise). 

Let A be the area of the pickup loop in Fig. ^, and let m be the unit 
vector normal to the loop. We let / and m be the real coefficients of the 
loop current I{t) = /e"^* and loop magnetic moment m{t) = me*'^*, such 
that m = Aim. We wish to compute the spectral density rh • Sb ' ''^ of 
the magnetic field traversing the loop. We reason as follows:^ 

m-SB''fn = Si^/A^, magnetic flux <^ = Am • B; 

= Sv/{<^A)'^, Faraday's law of Induction; 



= Ke{Z)/{ujAf 2kBT, thermal voltage noise per Mz. 

= P/{ujIAf 2kBT, because P = Re{Z)P; 

= m • r • Tn/{IA)'^ 2kBT, definition of T in terms of P; 

= rh - T • rh 2A;bT, substitute m = I Am. 

Since the result holds for arbitrary loop orientation m, and both Sb and V 
are symmetric matrices, it must be that Sb{^) = '^{^) ^^k^T, QED. 



A.l. Spatial correlations and gradients. This reasoning is readily ex- 
tended to encompass spatial correlations and gradients. Let n magnetic mo- 
ments [m] = {mi, m2, . . . , m„} be located at coordinates {iCi, £C2, • • • , x^}. 
Assuming a time dependence oc e*"^*, the dissipated power is calculable from 
the geometry and conductivity of the system — at least in principle — and will 
be of the general form P = [m] • [T] • [m] u^, where [T] is the obvious n-point 



''We write 2kBT in place of fiw coth(ftcj/(2feBr)) purely for economy of notation. 
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generalization of T. Then the n-point generahzation of ( |2.4| ) is simply 

(A.l) S[B] = [r]hu;coth(^^^ 

In this fashion all questions involving n-point spectral densities of magnetic 
noise can be reduced to the calculation of [F]. Although we will not do so in 
this article, the spectral density of magnetic gradients can be determined by 
considering the limiting case in which two coordinates approach each other. 

A. 2. Classical and quantum calculations. Students in particular may 
be troubled that our calculation of the dissipation tensor r(cj) will be en- 
tirely classical, via Maxwell's equations. After all, doesn't the title of this ar- 
ticle promise "the classical and quantum theory of thermal magnetic noise" ? 
The quantum part of the promise will be fullfilled in Appendix |^, to follow, 
in which we present a rigorous quantum theory of a spin coupled to a ther- 
mal reservoir, and show that the classical dissipative kernel determines the 
quantum dissipation also. 

A. 3. Solving the Maxwell equations. To calculate power dissipation 
within the slab, and thus determine r(a;), we seek solutions of Maxwell's 
equations in their standard form for time dependence oc e' 

(A.2) V X £; = -iuB V -D = pc 

V X H = j + iojD V -B = 

together with the constitutive and divergence equations for the current 

(A.3) j = {aE + f) V • j = -iujp^ 

Here B is the magnetic field, H = B/fi is the magnetic intensity, p, is 
the magnetic permeability, E is the electric field, D = eE is the electric 
displacement, j is the total current, e is the electric permittivity, pc is the 
charge density, a is the conductivity, and is the externally- applied dipole 
source current. 

A. 3.1. The scalar Maxwell equations. We begin our solution of the Maxwell 
equations with the help of a simplifying ansatz: that the irrotational part of 
E and D vanishes, i.e., the charge density pc is zero within the slabs and at 
their interfaces. Furthermore, we will assume that the displacement current 
iujD = iooeE is negligibly small compared to the conduction current aE. 
This approximation is quite accurate for all the cases we will consider, e.g., 
for room temperature copper at one GHz we have ue/a ~ 10~^]^ 

It is natural to write a divergence-free E field as a curl; we write E as 

(A.4a) E = -iujV X {ipn) 

^We remark that the irrotational ansatz and the neglect of displacement currents must 
be embraced simultaneously, because individually the ansatz and the approximation yield 
a mathematically inconsistent set of Maxwell equations. See, e.g., fell for a discussion. 



2ih 
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where V is a scalar vorticity, and we recall that h is the (constant) unit 
vector normal to the slab surfaces. Per our ansatz, ( A.4a| ) guarantees that 



V • -D = everywhere within a slab (boundary conditions will be considered 
shortly). In terms of ip, the Maxwell equation V X -E = —iojB becomes 

(A.4b) B = n (VV) - V(n • V^)) 

which, furthermore, guarantees V • -B = 0. The sole remaining Maxwell 
equation V X = aE + is then equivalent to the scalar wave equation 

(A.4c) V^V = WV' + / 

where is a scalar dipole source whose nature will be specificed shortly. 

It remains only to enforce the usual boundary conditions of electrodynam- 
ics: continuity of the normal components of B and D and the transverse 
components of E and H. These are equivalent to the following boundary 
conditions on ifj at the juncture of two slabs i and j 

(A.4d) = V'i 

(A.4e) (n • VV'i)//ii = (n • 

The problem is now reduced to the solution of the scalar wave equation 
(|A.4c ), subject to the between-slab boundary conditions (|A.4d ). The E and 



B fields are given by ( A.4a ) and ( A. 4b ), which now are purely constitutive 
equations. 

A. 3. 2. The scalar dipole source. For a (vector) magnetic dipole m located 
within a region of zero conductivity, it is convenient to specify the scalar 
source implicitly by j'^ = V'^tp^, with the source field tp^ given by 

47r V kc \x\+x-nj 



Inserting tp^ in the constitutive equation ( A.4b| ) yields for B 



(A.6) B(x) = ^ (3x0x-I)-m 

47r\x[^ 

which (by construction) is the dipole B-field generated by a magnetic mo- 
ment m. This justifies our specification of as the source field. The pos- 
sibility of writing a general dipole field in this form is the key step leading 
to the success of our scalar ansatz. 

As a technical point, the vector potential ip^{x)n describes an E-field 
which is singular along a line extending from a; = to infinity in the — n 
direction (away from the slab). As we will show in Section [A. 3. 6 , the singular 



portion of the E-field is irrotational, and it can be written as the gradient 
of a harmonic potential which disappears when we introduce a finite charge 
density on the surface of the slab. The dissipated power within the slab 
is not appreciably affected by the currents required to create this surface 
charge, so for the time being the singularity can simply be ignored. 
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A. 3. 3. Solving the scalar Maxwell equations. The scalar equation ( A.4c ) for 
ip is most easily solved in cylindrical coordinates {r,(p,z}. It is further 
convenient to work with the Bessel transform tpm{p^ z) defined bypl 



1 



2n roo 

(A.7) ^^{p^z) = ^l dip rdr e~'"''^J^{pr)i;{r,p,z) 

Jo 



27r 

whose inverse is 



m=oo 

(A.8) Hr,p,z)= ^ e^'"^ / dp p Jn,{pr)i^m{p, z) 

m=-oo •'O 

Here Jm is a Bessel function of integer argument. The scalar wave equation 
for '4>m{p,z) is 

(A.9) Q^^rnip, z) = zujfiaip„,{p, z) + f^{p, z) 

whose general homogenous solution is 

(A.IO) i^mip, z) = iPMe"' + i^Me^'' 



where we have introduced the complex wavenumber k = ^ iiofia + p^. 

Our next task is to explicitly solve for the Bessel coefficients il)^{p). We 
divide space into three regions {1,11,111} , where Region I is the nonconductive 
region that contains the dipole source, Region II is the conductive slab of 
thickness t beginning at a distance d from the source, and Region III is the 
nonconducting region on the opposite side of the slab (see Fig. |^). 

Within Regions I and III we assume vanishing conductivity o" = and a 
purely real permeability /Uq- In contrast, within Region II we allow both a 
and ji to assume complex values. Thus energy dissipation can occur only 
within Region II. 

Requiring that tpmip, — > as \z\ — > oo leads to the general solution: 
(A.ll) 

(P)e''^ + V'm(P> z) Region I (source); 

I z=d 

i>mip,z) = { v5m"^(/o)e'''^ + i^^- {p)e-^'' +z Region II (slab); 

i z=d+t 

Wm~ {p)^~''''^ Region III (far side). 



The Bessel functions Jm satisfy an ortliogonality relation 

rdr Jm{pr)Jm'{p'r) = S^m' Kp~ P)Ip 
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For illustrations of Regions I, II, and III, see Figs. |T|-|2lThe source term 
ip^{p, z) in Region I is found by substituting ( [A.5|) in ([A.?]) : 



(A.12) 
(A.13) 



/^O -p\z\ - 



Svrp 



e {x ± iy) • m 



where we have introduced Cartesian basis vectors x X y = h. Note the 
singularity at z = 0; this reflects our assumption of a point dipole source. 



Section A. 3.6 derives explicitly nonsingular expressions for B and E origi- 
nating from finite-sized current sources, however the present dipole result is 
adequate for our main purpose of estimating power dissipation in the slab. 



With the source term now specified, the boundary equations (A.4d) for 



the Bessel coefficients {V'o ; ' V'o '^o } '^^^ readily solved. We 
obtain for the m = coefficients: 



(A.14) 









_ Po 




47r 







g-2pd(g2H_i)(^2^2_^2 

_2e-{^'+'=)'='(Kp - k)Kp 
2e^^^-(p-'')'^{Kp + k)Kp 
4e(p+^)kKp 

n • m 



pUK^p"^ + A;2)(e2fct - 1) + 2Kpk{e'^''t + 1)) 

where K = p/po is the relative permeability of the slab. The m = ±1 
coefficients are directly proportional to the m = coefficients 



(A.15) 



{x lb iy) • m 
2(n • m) 



1+' 



When we calculate magnetic spectral densities, this simple proportionality 
will ensure that F oc ( i" + n ® n) at all frequencies and length scales. 

A. 3. 4. Energy dissipation. Our next task is to calculate the power P dissi- 
pated within Region II. From classical electrodynamics we have 



(A.16) 
P -- 



dV (j-E + H 



dB 



,v \ dt 
which with the help of the Bessel orthogonality relations becomes 
(A.17) 

°° roo rd+t / T™/',.^\ 

\i^m{p,z)\^ 



P = 27rV I pdp I dz ( p2t^2j^e(cr) + 



Im(^) \ 
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Carrying through the z-integration, we find the scalar magnetic dissipation 
coefficient ^{u>), as implicitly defined by P = lo'^T{uj) (m • (/ + n ® n) • m) : 

(A.18) r(a;) = ^ 1^ ( p^Re{a) + ^1^] iKpg-Sprf ^ 
4vr Jo V ) 



Re 



,2fct 



2fct 



\{K'^p^ + A;2)(e2'=* - 1) + 2Kpk{e^^^ + 1) 



where we recall that k = \J iu>fia + p'^ and K = fi/ fiQ. In simplifying this 
result we have kept only the leading and next-to-leading terms in Im(/i), 
i.e., we have assumed K is predominantly real. 

The case of thermal magnetic fluctuations measured at a point midway 
between two identical conducting slabs can be solved by exactly similar 
methods. It is only necessary to replace the Region I homogenous term 
ipl^{p)e''^ in ([A. 11 ) with the symmetrized term ip}^ (p) cosh.{pz) . The re- 
sulting magnetic dissipation coefficient is 



(A.19) r' 



p^Im(;u) \ 



■2pd 



2Re [iA;*(e2'=* - l)(e2^** + l){K^p^ - k'^)/{ujpa) + pKje^'^* - Ip] 
|(/t2 + K^pi + e-2p'i(fe2 - K2^2))(e2fct _ 1) + 2Kpk{e^''t + l)f 

Note that the double-slab integrand is just twice the single-slab integrand, 
slightly modified by an additional term tx e~^'"^ in the denominator. For the 
special case K = 1 (the case of primary interest) this additional term can be 
shown to be non-leading for all values of {p,d,t,a,uj}. A numerical survey 
confirms that r'(ci;) ~ 2r(ti;) to within very good accuracy, as described in 
Section 2.1. 



A. 3. 5. Asymptotic limits for paramagnetic and diamagnetic slabs. The three 
asymptotic expressions in (|2.5|) are obtained by evaluating ( A.18 ) for K = 1 
and Im(/i) = in the following limits: 



' (1 - e-2p*) 



(A.20) r(tj)~Re((j)^ rdpx<> 

4vr Jo 



8g2pd ' 
^3^3g-2pd 

2cos(7r/4-(/./2)' 
p^X'^it - 4pA2 sin<? 



\a\ — > oo; 



-, A;t < 1 «; pd. 



The integrations are readily carried out, with the results as in (p.5|). 



A. 3. 6. The irrotational electric field. Now we have only one chore remain- 
ing: repairing the irrotational singularity in the E-field mentioned in Sec- 
tion A. 3. 2 , and verifying that the repair does not alter the energy dissipation 
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just calculated. In the end, we will obtain expressions for the E- and B-fields 
which are closed- form and explicitly finite. 

In outline, the repair is accomplished by adding a purely irrotational 



E-field to Region I — thereby relaxing the ansatz of Section A. 3.1 that E 
is purely rotational.[] With the added E-field, Faraday's law ^ X E = 
—icuB is still satisfied, because the added E-field is irrotational. Because 
the added field is confined to Region I, which is nonconductive, and because 
displacement currents iioD are negligibly small, Ampere's law V X -ff = 
aE + j'^ + iujD ~ is still satisfied. Gauss's law V • = is still satisfied 
because the added E-field is the gradient of a harmonic potential. Physically, 
the added field is generated by surface charge on the slab, and the within- 
slab currents necessary to sustain the surface charge are negligibly small 
compared to the eddy currents already computed. Thus, the added E-field 
does not alter the magnetic dissipation coefficient r(a;) already computed. 
The remainder of this section is devoted to proving these assertions. 

We begin by noting that the E-field singularity in Region I physically cor- 
responds to the field induced by an end-to-end line of electric dipoles along 
the — n axis. It is easy to verify that such a singularity can be cancelled by 
the addition of an (irrotational) E-field which is the gradient of the harmonic 
potential (^{{x), 

fA ^1^ ^ . no x-{hx m) 

(A. 21) 'Pi\X) = ■ 



47r -|- a; • n 

To satisfy the boundary conditions between Regions I and II, we must also 
add the E-field of the image potential 

(A.22) 4>l{x) = -(f>\{2dh - x) 

in order that E^ = V(0| + (/'2) normal to the slab surface, thus ensuring 
that the transverse component of E remains continuous at the boundary. 

To sustain the E^, a charge density = efi- E^ must exist on the surface 
of the slab, and this surface charge must in turn be created by a normal-to- 
surface current within the slab = iujp^h. Relative to the intraslab eddy 
currents current aE which are our main concern, it is easy to show that 

is of order uje/a ^ 10~^ (see the discussion following ( [A.2| )) and thus is 
negligible. 

Putting all these pieces together, we are now ready to specify the fields 
in Region I in a manifestly finite manner, assuming only that the spatial 
size of the current source is reasonably small compared to the distance d 
separating the source from the slab (like the current loop in Fig. ^)J^ 



We recall Helmholtz's theorem: that any finite vector field that vanishes at infinity 
can be written uniquely as the sum of an irrotational part (the divergence of a scalar 
potential) and a rotational part (the curl of a vector potential) 

**Also, we continue to ignore the displacement current iujD in Ampere's law V X = 
j + iuD\ this is equivalent to the near-field assumption ud <Si c, where c = (/ioeo)^^'''^ is 
the speed of light. This near-field approximation is very well satisfied for most practical 
problems in quantum-coherent engineering. 
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From the localized current source, we compute a vector potential A{x) 
in the usual manner 

(A.23) A(.) = f /rfV^i^-^^ 

^ Att Jy \x - x'\ 47r |a;|2 

where j{x') is the (finite) current density within the (finite) source. From 
the far-field limit of A{x) we determine the dipole moment m of the source, 
and under the assumption that the size of the source is small compared to d, 
we use m to compute the back-reaction vorticity field '^^(a;) via equations 
( A.$ - A.14 ) and the back-reaction irrotational field via ( A. 22 ). With the 



help of the identity —iu}A{x) = —ico'V X {tp^^x)) + ^(j)\{x), which is exact 
for a point dipole current source, the net B-field and E-field in Region I can 
be written as follows 

(A.24a) B{x) = V X A{x) -V{h-V7p\x)) 

(A. 24b) E{x) = -iLoA{x) - iwV X {i^\x)h) + V(pl{x) 

where A{x), ip^{x), and (p^ix) are individually nonsingular. 

Physically speaking, A{x) describes the E- and B-fields that are generated 
directly by the source current, ip^{x) describes the E- and B-fields generated 
by the induced currents within the slab (treating the source current as a 
dipole), and </'2(^e) describes the E-field generated by surface charges on the 
slab (also treating the source current as a dipole). These expressions were 
used to generate the field geometry in Fig. 0. 



Appendix B. Quantum decoherence 

In this section we present an exactly solvable microscopic model of a 
spin-^ particle magnetically coupled to a heat bath. Our goal is derive the 
expressions for Ti, Ta, and Tip quoted in (^)-(^). 

We follow Ford, Lewis, and O'Connell in modeling the heat bath as a 
collection of independent harmonic oscillators the reader is referred to 
this article for a thermodynamic justification of this model. We will adjust 
the bath parameters so as to reproduce the known spectral density of mag- 
netic fluctuations, which we computed in the previous section. The quantum 
Langevin and optical Bloch equations — which yield closed-form expressions 
for Ti, T2, and Tip — will then be uniquely determined. 



B.l. The thermal reservoir Hamiltonian. The Hamiltonian H of the 
spin/heat bath system is taken to be 

(B.l) H = -7B0 -s + Y^ [p] + {qj - (3jhj • sf) 

j 

Here Bq is the polarizing field, 7 is the gyromagnetic ratio of the spin, and 
s is the angular momentum operator of the spin, satisfying commutation 
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relations [si,Sj] = ihcij^Sk- The heat bath is described in terms of oscilla- 
tors with resonant frequency Wi, whose dynamical coodinates {pi, qi} satisfy 
[Qi^Pj] = ifi^ij, which are coupled to the spin with strength (3j. 

The physical nature of the heat bath variables need not be otherwise 
specified. Reader may optionally conceive them as, e.g., thermal excitations 
of conduction band electrons, magnons in a ferromagnet, excited states of 
the vacuum in field theory, or in general as excitations of whatever heat bath 
is furnished by a given quantum technology. 

We begin by solving the equation of motion of the heat bath variables. 
In the Heisenberg picture we have|^ 

(B.2) qj + Lu'jqj = PjUj'jhj • s 

Following Ford et al, we write the formal solution to these equations as 
(B.3) qj{t) = q^{t) + Pjfij • s{t) - f dt' cos(wj(t - t')) PjUj • s{t') 

J — oo 

Here q^{t) is a homogenous solution to the equations of motion (correspond- 
ing physically to the evolution of the heat bath in the absence of the back- 
action of the spin). With the help of this result, the spin operator equation 
of motion s = [s, H]/{ih) can be readily cast into the form of the following 
quantum Langevin equatior^^ 

(B.4) sit) = -7S0 X s{t) - 

^-f{B{t) X s{t) - s{t) X B{t)) + 

\ [{C.s{t)) xs{t)-s{t) X {C-s{t))] + 

/ dt' - [{G{t - t') ■ s{t')) X s{t) - s{t) X {Git - t') ■ s{t'))] 
J -00 ^ 

Here -B(t) is a fluctuating thermal magnetic field which depends only on the 
homogenous heat bath operators (t) 

(B.5) S(i) = J^T^Wi^j^W^i 

i 



^We recall that the Heisenberg equation of motion for an arbitrary operator 0(i) is 
0(€) = ^ ,H\I (iK). The results of this section all follow directly from this equation, 
using only the commutation relations for {qi,pi, s}, plus the fact that commutators in the 
Heisenberg representation are the same as those in the Schroedinger representation. 

^"As a help to students, we note that the quantum Langevin equation can be written 
in several equivalent forms, which arise because the commutation relation [qj{t), s{t)] = 
allows the heat bath interaction to be written in any of the following fully equivalent forms 

q,{t)sit) = s{t)q,{t) = + s{t)q,{t)) 



Upon substituting (B.3) for qj{t), these various forms yield fully equivalent (but super- 
ficially quite different) Langevin equations. We have chosen to work with the fully sym- 
metrized Langevin equation; none of our main results depend on this choice. 
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and G{t — t') IS a dissipative kernel 

(B 61 G(t -t') = i "^^^^ ~ '^^ * ^ 

^ ' ' ^ ^ [0 t<t' 

Note that G{t — t') is a pure c-number matrix which is independent of the 
temperature of the heat bath. It is the matrix generahzation of what Ford 
et al. call the "memory function" of the heat bath. 

The anisotropy tensor C is a symmetric c-number tensor given by 

(B.7) ^ = XI ^^^J '^i ® ""J 

3 

Physically speaking, C describes the interaction of the particle with image 
currents in the nearby slab. This reflects a very general and physically real- 
istic property of the independent oscillator model: the dynamical equations 
of a particle are renormalized by its interaction with the thermal reservoir. 

For the special case of spin-| particles C has no dynamical effects, due to 
an identity satisfied by spin-^ operators (and not by higher-spin particles) 



2 

where I is the identity operator. This identify ensures that the C-dependent 



(B.8) SiSj + SjSi = — I 



terms in (B^) vanish identically; we emphasize that this occurs only for 
spin-i particles. 

So far, we have not used the fact that the heat bath variables are in 
thermal equilibrium; we now use this fact to compute the spectral density 
Sb of B{t). Not surprisingly, we will discover that Sb has a functional 
form that is closely related to the dissipative kernel G. 

The calculation is straightforward. Denoting a thermal ensemble aver- 
age by (...)*, we follow Ford et al. in noting that the correlation of the 
homogenous heat bath operators {qf} is of the simple form 

(B.9) {qHt)q^{t'))t = \ hcoih cos(a;,(t - t')) 

Combining this with ( B.5D yields an explicit expression for SBi^) 



(B.IO) 

SB{io)= / dTe-'''^-{B{t)^B{t + T) + B{t + T)^B{t))t 

fij (g) fij hcoth. ( f^^ ^ [ dt' e"*"^^ cos(a;or) 
r \2kBTJ 



-oo 

3 
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2^ ^3 ® ^3 ^ COth J [H^ - ^j) + ^{^ + ^j)) 
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If we compare this result with the real part of the Fourier transform of the 
dissipative kernel 

/oo 
dt'e~'^^^-^'^G{t-t') 
-oo 

= ^ ^ LOjPj flj (g) flj [5{uj - UJj) + 5{ijJ + UJj)) 
j 

it is apparent that 

(B.12) Sb{uj) = -f-^Re G{lo) hw coth 

In terms of the magnetic dissipation matrix T{uj) computed in the previous 
section, this takes an even simpler form 

(B.13) -i'^T{uj) =ReG{uj) 

This is the fluctuation-dissipation theorem as it applies to particles of arbi- 
trary spin. For our purposes an alternative relation is even more useful: if 
r(a;) is known, then the dissipative kernel G{t — t') is given explicitly by 

„2 /-oo 

rf^eMt-t')r(a;), t>t'; 



(B.14) G{t-t') 




t < t' 



as may be verified by writing T{uj) in terms of Sb{^^) using (2^), then 
writing Ssi^) in terms of heat bath variables using ( plop , then carrying 
out the integration and comparing with ( |B.6| ).F^ 

It may occur to the reader that by virtue of ( |B.13| ) and (B.14), knowledge of Re G(lj) 
suffices to determine Im G{lo). This insight leads directly to a Kramers-Kroenig relation: 

/oc 
-OO 

/■OO roc ^ 2 

= lm dr duj' e''-'^ ~"^^-^r{uj') 

Jo i-oo 
-j^ poo roo 

= — dr dcj' sin ((lu' — Lu)r) KeGiuj') 

Jo J-oo 

which is the Kramers-Kroenig relation for particles of arbitrary spin interacting with 
thermal magnetic noise. Most textbooks present such relations in a more compact but 
less obviously finite form, which can be obtained by rewriting the above result as 

1 dr.' U ' ^ ^ Re(G(a;') - GH) 



— I dr —— I duj cos ((cu — a;)r) 
Jo dTj_^ 



where the G{uj) term is introduced to regulate the singularity at uj' = ui, and the co' inte- 
gration contour is adjusted in the complex plane so that the net contribution of this added 
term is zero. The r integration then becomes trivial and the remaining u' integration is 
the Kramers-Kroenig relation in its traditional form 

(B.15) imGic.) = - Pf^dc''^^^^ 
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B.2. The master equation. We now direct our attention to solving the 
quantum Langevin equation ( [B.4D . As it stands, this equation is too comph- 
cated to solve directly]^ In particular, B{t) is an operator-valued function 
of heat bath variables whose detailed dynamical behavior — describing every 
minute fluctuation within the thermal reservoir — we neither know nor wish 
to know. To make progress, we will average over a thermodynamic ensemble 
of heat bath variables, and study only a coarse-grained time derivative of 
s{t). Simplified equations obtained by this general strategy are called mas- 
ter equations — our goal is to derive a master equation from the quantum 
Langevin equation. 

We begin by defining a rotation matrix R{t), satisfying R{t) = 'jBq X 
R{t), in terms of which we define rotating- frame quantities 

s,ot{t) = R{t) ■ s{t) 
B,ot{t) = R{t) • B{t) 
G,ot{t-t') = R{t)-G{t-t')-R^{t) 

It follows that s{t) = R^t) • (srot(i) - jBq X Srot(i))- 

Neglecting 0(s^q^-) terms — this is the rotating- frame approximation — and 
specializing to spin-i particles so that the anisotropy tensor C does not 
enter, the quantum Langevin equation ( [B.4| ) becomes 

(B.16) S,ot{t) = -^7(Srot(t) X Srotit) - Srot(0 X Brot{t)) + 

^ J dt' ((Grot(t - t') ■ (7B0 X Srot(t'))) X Srot{t') - 

Srot(t') X {Grotit - t') ■ (7S0 X Srot(t')))) 

Recalling that Gmt{T) is the temperature-independent "memory function" 
of the heat bath, we now assume it decorrelates rapidly compared to the rate 



where P is a principal value. 

^■^ Students (and engineers) may wish to reflect that the quantum Langevin equa- 
tion (B.4) is easy to solve in principle — it can be numerically integrated quite readily 
(the integration is quite straightforward to program in languages like Mathematica; this 
is a good exercise for a student) . Heisenberg operators are be represented numerically by 
matrices of complex numbers, with operator addition and multiplication achieved via ordi- 
nary matrix addition and multiplication. As time goes on, the quantum Langevin equation 
increaingly entangles the spin operator s{t) with the thermal reservoir basis states — thus 
leading, as expected, to fluctuation, dissipation, and entanglement of the spin operators. 

The only practical difficulty with this program is that a reservoir of n particles, each 
having m quantum states, requires m" quantum basis states, with each Heisenberg oper- 
ator stored as a time-dependent m" x m" Hermitian matrix. For n as small as two and 
m as small as twenty — too small to describe a realistic thermal reservoir — the storage and 
multiplication of these exponentially large matrices is enough to overwhelm any classi- 
cal computer. Quantum computers were invented in part to overcome this fundamental 
quantum simulation problem. 



FLUCTUATION, DISSIPATION, AND ENTANGLEMENT 



30 



at which Sj.ot{t) varies. We can then ignore the time-dependence of s,.ot{t) 
inside of integrals, and make the fohowing substitution 

(B.17) r dt' G,ot{t - t') = 7'rrot(0) 

J —oo 

Here the rotating-frame magnetic dissipation matrix FrotC'^) is related to 
the laboratory- frame matrix T{uj) via 

(B.18) rrot(w) = [6 ® S] tr [& 6 • T{lo)] + 

[l-b^b] itr[(/-6®6) •r(wo)] 

with 7S0 = iOob. Similarly, Ss^oti^) related to Ssi^) via ( 1 2. 8| ). Then, 
with the help of some straightforward algebraic manipulations ,p[the Lange- 
vin equation for spin-^ particles can be cast into the form of a master equa- 
tion, which turns out to be the following Block equation 

(B.19) ^ (srot(t)) = - ^ 6 ® • ((srot(t)) - {so)) 

iI-b^b)-{s,ot{t)) 



Here Ti and T2 are as given in ( |2.9| ), and we have wrapped the spin operator 
in an ensemble average (...) in order to express the Bloch equations in 
terms of the (c- number) vector spin polarization {sj-^t), as is traditional. 
The equilibrium spin polarization (sq) is found to be 

(B.20) ^ I .a„>. (gl) 6 

which agrees with the well-known thermodynamic expression for spin-^ po- 
larization; this provides a consistency check of the Langevin formalism. 

The calculation of the spin-locked relaxation time Tip is similar, differing 
only in that transformation to a doubly rotating frame is required. The 
result is as given in (p.9D; details of the derivation will not be given. 



^■^The derivation involves integrating (B.16) to first order in Grot and second order 
in Srot, then substituting (B.17). The rotating-frame commutation relations are then 
invoked to simplify the cross-terms: 

[(Srot(i))i, (Srot(t))j] = ifieijk{Srot{t))k, 

These are valid because spin commutators are form-invariant under transformation to the 
Heisenberg picture and to the rotating frame. Further specializing to spin-i implies 

where I is the identity operator; these identities also are form-invariant. Finally, when we 
average over a thermal ensemble, we have 



dt' i {Brot{t) <E) B,ot{t')+t^ t') = Ss„,(0) 



as well as (trivially) 
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We close by noting that our spin relaxation rates differ from the oscillator 
relaxation rates of Ford et al. Q in one physically important respect: Ti, T2, 
and Tip are strongly temperature-dependent, while Ford et al. predict an os- 
cillator quality Q that is independent of temperature. Formally, the reason 
for this difference can be traced to the Langevin equation ( B.16| ), in which 



the dissipative kernel Grot is temperature-independent (this is true for both 
spins and oscillators). However, when this Langevin equation is integrated 
to second order to obtain the Bloch equations, the temperature-dependent 
fluctuations B{t) also contribute to spin relaxation via the commutation 
relation [sj, Sj\ = iheij^Sk- In contrast, the oscillator operators p and q have 
a pure c-number commutator [q,p] = ih, and in consequence the fluctuat- 
ing force exerted by the reservoir does not directly contribute to oscillator 
relaxation. 

The difference between spin relaxation and oscillator relaxation can also 
be understood physically. If we imagine a classical spin that is subject to 
a randomly fluctuating magnetic field, and no other physical influence, the 
average spin orientation (s) necessarily will relax toward (s) = 0. So to the 
extent that a higher-temperature reservoir creates stronger magnetic fluc- 
tuations, spins will relax more rapidly at higher temperature. In contrast, 
a randomly fluctuating force does not change the average trajectory of an 
oscillator at all — provided the average force is zero — because oscillator equa- 
tions of motion are linear, unlike spin equations. So to the extent that a 
thermal reservoir has a temperature-independent dissipative kernel — like the 
independent oscillator model — the oscillator quality Q will be independent 
of temperature. 

Appendix C. Fluctuation-dissipation-entanglement theorems 

Now we will prove fluctuation-dissipation-entanglement theorems for both 
spin-^ particles and harmonic oscillators in contact with a thermal reservoir. 

C.l. A spin-i fluctuation-dissipation-entanglement theorem. Sec- 



tion |2.3| has already provided the definitions we need to prove the fluctuation- 
dissipation-entanglement theorem, and Appendix |b| has carried through 
many of the needed calculations. We need only organize our reasoning as 
follows. 

Specializing to spin-^ particles, we write the total Hamiltonian ( [B.l[) in 
the usual form H = Hq -\- V, with the perturbing Hamiltonian V given by 

(C.l) V = -Y,^jPjQj'n'j-^ 

j 

The exact ground state |^'o) of H can be calculated order- by-order in V as 



in the notation of Landau and Lifschitz |1C]. The formal expression for 



I^'q^^) is quite lengthy, but fortunately we have £|\I'q''^) = by construction 
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( Un\ ), so that cross terms like (^'q^^ |£|^'q'''') vanish. In consequence, the 
leading contribution to the entanglement comes entirely from (^'q^'*|£|\I'q^^), 
for which we have the well-known expression 



(C.2) 



En — Eri 



The fluctuation-dissipation-entanglement theorem ( 2.12 ) asserted in Sec- 
tion |2.3| can now be readily derived 



(C.3a) E{p,oj^) = (vI/o|£|vI/o) ~ (^ri£|^|,' 



(C.3b) 
(C.3c) 

(C.3d) 
(C.3e) 



8(u;o + '^'jY 



tr [{I-p®p)- {hj ® )] by Q 



An. 



6{uj — iOj) 



tr [(/ -p®p) • {hj (g) hj)] 

^ 02{u;tr[(I-p(g)p).Re(G(u;))];u;o} 
02{c^tr[(/-p(g)p) •r(cj)];a;o} 



h 

in 



by ( iBlll ) 
by (iH) 



where ([03^) was obtained by introducing tj as a dummy variable of inte- 
gration over 6{lu — ujj)- 

C.2. An oscillator fiuctuation-dissipation-entanglement theorem. 

A similar fluctuation-dissipation-entanglement theorem exists for harmonic 
oscillators. We specify the Hamiltonian H of the oscillator/thermal reservoir 

as 

(C.4) H = ^uo {p" + g') + E l^'^f^i + - 

j 

Here q and p are generalized oscillator coordinates satisfying [q,p] = ih, 
whose physical nature is not otherwise specified. As shown by Ford, Lewis, 
and O' Council Q, the Hamiltonian ( p.4| ) implies the following quantum 
Langevin equation in the Heisenberg picture 

ft 



(C.5) 



q{t) + cjgg(t) +uJo[ dt' iJL{t - t')q{t' 



fit) 



where f{t) is a fluctuating thermal force, and ^{t) is a dissipative "memory 
function" given by 



(C.6) 



^^{r) 



Yl,j ^j0j cos{ujjT) for r > 0; 
for T < 0. 
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The Fourier transform of /x(t) 

/oo 
-oo 

has a real part given exphcitly by 

(C.8) Re[/2(a;)] = \Y^ ^il^l (^^^ " + ^^"^ + '^i)) 

3 

which is related to the spectral density Sf{uj) of the thermal force f{t) by 
(C.9) Sf{u;) = Lol Re[A(a;)] hio coth 

This is, of course, the celebrated fluctuation-dissipation theorem for har- 
monic oscillators. 

Now we have all the pieces we need to calculate the zero-temperature en- 
tanglement from the dissipative kernel jj,{uj). By reasoning exactly analogous 
to the spin-^ case we find 

(C.lOa) (^'o|e|*o)=^K'^|£|^[,'^> 



(C.lOb) =^ 



J 

oo 



4(u;o + 



(C.lOc) = / dio- 5:^^—7^ Re[/2(u;)l 

Jo ('^o + 

(C.lOd) =^ (J52{^Re[/i(c^)];c^o} 

This is the general form of the fluctuation-dissipation-entanglement theorem 
for harmonic oscillators. 

As in Section 2.3. Ij for spin-^ particles, oscillator entanglement can be 
evaluated approximately. We stipulate that jl{uj) ~ l/Q for w G {0,uJc), 
with Q the oscillator quality and lOq the thermal reservoir's cutoff frequency 

3> wq. Then we find to leading order in loq/uJc 

(C.ll) (*„|e|*„, . 

This approximate result is not new. Li, Ford, and O'Connell |^, eq. 14] 
have calculated a zero-temperature energy shift that directly implies ( C.llD , 
in the context of their reply to a critique by Senitzky [^] of a still-earlier 
analysis of energy balance in dissipative systems |15]. However, they did not 
interpret the energy shift as an entanglement relation, and the generality 
and rigor of the link between dissipation and entanglement via Stieltjies 
transforms was not pointed out. 



FLUCTUATION, DISSIPATION, AND ENTANGLEMENT 



34 



C.3. Oscillator renormalization. For spin-^ particles, we have seen that 
the thermal reservoir interactions renormalize the dynamical equations of 



the field operators, via the anisotropy tensor C (B.4,B.7). We will now show 
that a similar renormalization occurs for oscillators. 

To begin, we remark that in quantum-coherent engineering, renormal- 
ization is more than an abstract concept. At least in the context of mag- 
netic resonance force microscopy, we will see that renormalization effects are 
large — they play a key role in experimental protocols — and are of consider- 
able practical consequence. 

Next, we note that oscillator renormalization is much easier to analyze if 
we recognize at the outset that the operators q and p and the frequency loq 
which appear in the independent oscillator Hamiltonian ( p. 4]) are already 
renormalized. We know this from the quantum Langevin equation ( |C.5| ), in 
which ujQ appears as the resonant frequency of the oscillator after it has been 
renormalized by interactions with the thermal reservoir. Our task, therefore, 
is to deduce an unrenormalized "bare" Hamiltonian from the renormalized 
"dressed" Hamiltonian ( p.4| ). A great virtue of the independent oscillator 
model is that this calculation can be carried through quite easily. 

We define an unrenormalized frequency uj'q by 

(C.12a) u;'o^[u;o{u;o + Ej^jP])f^ 

and unrenormalized operators q' and p' by 

(C.12b) q'^qiuyuo)'/^ 
(C.12c) p'^p(^oM)'/' 

and unrenormalized reservoir couplings by 
(C.12d) ^ (u;oM)i/2 

By construction, the unrenormalized operators q' and p' satisfy the canonical 
commutation relation = ih. Furthermore, the unrenormalized param- 

eters Wg and are such that the Hamiltonian — when written in terms 

of unrenormalized quantities — takes the desired "bare" form H = Hq + V , 
with Hq and V given by 



(C.13a) Hq = \j,{p'' + q'') + J] \u:,{p] + q^) 

3 

(C.lSb) V = -Y,^,q,f)'^q' 



The virtue of the "bare" Hamiltonian Hq is that it clearly shows what hap- 
pens when the oscillator-reservoir couplings {P'j} are turned off: the resulting 
"bare" oscillator dynamics are described in terms of unrenormalized canon- 
ical operators q' and p' and resonant frequency uj'q. 
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The physical meaning of — and the necessity for — renormahzation is made 
clear by the following real-world example. In magnetic resonance force mi- 
croscopy (MRFM) experiments, the cantilever resonant frequency is rou- 
tinely monitored as the tip of the cantilever is brought near to a sample. 
As the tip nears the sample, the Van der Waals potential (equivalent to 
the Casimir effect) increasingly acts to reduce the cantilever's spring con- 
stant, which renormalizes the resonant frequency ojq to a lower value than 
the "bare" value u>q, consistent with ( |C.12a ). Renormahzation effects in 
MRFM are so strong that is commonplace for the renormalized frequency 
to pass through zero and become imaginary, in which case the cantilever 
becomes unstable and the tip "snaps in" to the sample. Observation of the 
renormalized frequency provides a vital MRFM technique for monitoring 
the tip-sample separation without ever touching the sample. In our exper- 
iments, we frequently operate at frequency shifts in excess of 100 Hertz, so 
that renormahzation is by no means a small effect. 

Concomitantly, the cantilever's quantum ground state is also renormal- 
ized: far from the sample the ground state satisfies (g' + «p')lV'o) — Oj 
while near to the sample the ground state satisfies {q + ip)\ipo) = 0. Physi- 
cally speaking, the renormalized ground state |^o) and the unrenormalized 
ground state {i/jq) are squeezed relative to one another. Although deliber- 
ately squeezing cantilever states is not yet routine practice, it is not too 
far-fetched to imagine that someday it may be. 

We close this section by remarking that — as is typically the case in renor- 
mahzation theory — the inverse problem of computing the renormalized pa- 
rameters from the bare parameters is not solvable in closed form.|^ Instead, 
the renormalized parameters must be calculated perturbatively, order by or- 
der in {Pj}, and convergence is not guaranteed. In fact, perturbative renor- 
mahzation is guaranteed to fail when the bare Hamiltonian has a spectrum 
that is unbounded from below, because we are trying to convert the bare 
Hamiltonian into a manifestly positive-definite renormalized form. Such di- 
vergences are perfectly physical: they are realized experimentally whenever 
a cantilever tip snaps in. Thus, all aspects of renormahzation theory — 
even perturbative divergences — acquire practical significance in quantum- 
coherent engineering. 



C.4. Renormalization-dissipation relations. Now we show that renor- 
mahzation parameters can be calculated from the same dissipative kernels 
that control fiuctuation, dissipation, and entanglement. In our nomencla- 
ture, these relations do not qualify as "theorems" because the relationship 
is not invertible — the form of the dissipative kernel cannot be inferred from 
measurements of renormahzation. 



^■^To see this, students should try to invert (C.12a-C.12d) to express the renormalized 
variables {g,p, a;o,/3j} as closed-form functions of the bare variables {q' ,p' ,u)'o, Pj}. 
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For oscillators, renormalization is entirely specified by the ratio of the 
unrenormalized frequency ujq to the renormalized frequency loq, per ( p. 12a - 
C.12d). From ( |C.^ ) and ( C.12a| ) we express the frequency shift in terms of 
the dissipative kernel Re[/i(a;)] 



/2 



(C.14) 



TT 



1 + - / duj Re[fL{u;)] 



For particles of arbitrary spin, renormalization is entirely specified by the 



anisotropy tensor C, per the quantum Langevin equation (B.4). From 
amd ( |B.ll ) we find the anisotropy tensor in terms of the dissipative kernel 



(C.15) 



C 



-/ dujRe[G{u)] 



We see that renormalization effects are temperature-independent; this is 
consistent with a careful reading of the independent-oscillator literature [0. 
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Figure 1. Conductive slab geometry. For a realistic rendi- 
tion of B-field and E-field geometry, see Fig. |2[ Th is figure 
defines Regions I, II, and III as referenced in ( A.ll|) . 
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Figure 2. The B-field, E-field, and intra-slab currents — 
drawn to accurate scale — created by external excitation of 
a finite-sized current loop, as specified by ( |A.25 - A.24). The 
field geometry has a ready physical interpretation: the coil 
excitation current generates a B-field (solid ribbons), whose 
time-rate-of change induces a toroidal E-field (circumferen- 
tial dashed ribbons), that in turn induces currents within the 
slab (streamlines on the slab surface) which create a return 
flux loop (broad dashed ribbon) linking the induced slab cur- 
rents back to the coil. The return flux induces a phase-lagged 
coil voltage equivalent to that of a resistive impedance, thus 
ensuring that the energy externally supplied to excite the coil 
balances the energy dissipated in the slab. 
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Figure 3. Predicted spin relaxation rates in magnetic res- 
onance force microscopy and in quantum computing, as a 
function of conductivity. 



